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The conservation equation for electrons in a laminar or turbulent flow has 
been used to determine a criterion for breakdown. The results can be used 
to define a characteristic length, L 8 , and time, t s , which determine the affects 
of the flow on the power required to break down the gas. The theory has been 
compared to the experimental results of Buchsbaum and Cottingham in hy- 
drogen with reasonable agreement in the laminar region. In the presence of 
turbulence, the results will depend on the velocity dependence of the tur- 
bulent diffusion coefficient for the electrons. 

I. INTRODUCTION 

The physics governing the breakdown of a stationary gas (no flow) 
has been discussed at great length by a number of authors. For some 
recent studies the interested reader is referred to a text by Brown 1 and 
papers by Buchsbaum 2 and by McDonald. 3 The effects of gas flow on 
breakdown has been studied experimentally for a laminar flow of less 
than 100 meters per second by Skinner and Brady 5 and for the case of 
laminar and turbulent flows by Buchsbaum and Cottingham. 4 

Theoretically, the effects of gas flow can be studied by noting that 
the transport of electrons and ions will depend on the flow as well as the 
laminar or turbulent diffusion. If the flow is laminar and the electron 
density such that the diffusion is ambipolar, then the average drift 
velocity of the electrons and ions are equal and given by 

v = -D a (Vn/n) + V (1) 

where D„ is the ambipolar diffusion coefficient, n, the electron density 
and V, the flow velocity. The ambipolar diffusion coefficient may, in 
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turn, be related to the diffusion coefficient for the ions, D, , and the 
temperature of the ions, 7\ , and electrons, T e , by 

D.-D<{1 + Te/Ti) (2) 

provided the mobility of the electrons is much larger than that of the 
ions. 

If the flow is turbulent, then the ambipolar diffusion coefficient must 
be replaced by the turbulent diffusion coefficient, D T . 

The theory to follow will contain ratios of the form 

V/2D and V 2 /2D, 

where V is the average velocity in a turbulent flow or simply the velocity 
in a laminar flow. The first can be used to define an inverse length, L, , 
such that 

V/2D = L.~\ (3) 

and the second an inverse time, t s , such that 

V 2 /2D = tr 1 . (4) 

The characteristic time, t, , and length, L a , are a measure of the extent 
to which electrons are removed by the flow. These may be compared to 
the removal of electrons by ordinary diffusion to the walls with the 
characteristic time, t D , given by 

D/A 2 = t D -\ (5) 

where A is the size of the container. If 

t, » t D , 

then it is reasonable to expect that the effects of the flow on breakdown 
are negligible. If, in turn, 

t D » t e , 

it is reasonable to expect that the effect of ordinary diffusion to the walls 
is negligible as compared to "sweeping" as an electron removal process. 
In this case, breakdown will be controlled by the flow. 

II. GENERAL THEORY 

The effects of gas flow on breakdown can be considered theoretically 
by noting that the transport of electrons, on the average, is due to the 
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mean flow as well as the laminar or the turbulent diffusion. In the ab- 
sence of secondary emission from the walls, the appropriate conserva- 
tion equation for the electrons is given by 6 

?? + V- (nV) £* (vi - va)n + an 2 + V- (DVn) (6) 

at 

where 

n = the electron density in a laminar flow or the ensemble average 
of the electron density (at any given time) in a turbulent flow. 
V = the velocity in a laminar flow or the ensemble average of the 
gas velocity in a turbulent flow. (Experimental conditions 
usually justify a time average in this case.) 
Vi = the ionization frequency 
v a = the attachment frequency 
a = the recombination coefficient 

D = the diffusion coefficient for the electrons (laminar or turbu- 
lent). 
Consider the geometry of Fig. 1. The originally ambient gas is flowing 
between two parallel planes separated by a distance d. Two grids are 
placed normal to the direction of flow and separated by a distance I. 
For this geometry, (6) reduces to, 
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If it is further assumed that the distance I in the direction of flow is not 
excessive, such that, 

%* = for 0^2/ 2g J, 
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Fig. 1 — Cross section of flow geometry. 
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t8 = l[^] {1 + Te/Ti] 7~ (17) 

or by 

t -§[$]u + ivjw£. as) 

As expected, both L s and t. approach infinity as the flow velocity, V u , 
approaches zero. 

The magnitude of the characteristic sweeping length, L, , must be 
compared to the length of the discharge, I, in the direction of flow, y. 
If L, » I, then the term (V y /2D)y, may be neglected. 

The magnitude of the characteristic sweeping time, t, , must be com- 
pared to the characteristic time for electron loss by diffusion to the 
walls, t D ■ If t D » U , then diffusion to the walls may be neglected as an 
electron removal process. If t, » Id , then sweeping may be neglected. 
It is evident that since 

U ~ V~ 2 (19) 

then, t a » t D , for low flow velocities. As expected, sweeping may be 
neglected as an electron removal process at low velocities. 

2.2 Turbulent Flow 

When the flow of the gas is turbulent, it is only reasonable to expect 
that the diffusion coefficient, D, will increase provided the electron den- 
sity is sufficient to insure ambipolar or near ambipolar diffusion in the 
absence of turbulence. While this may be confusing, it is of importance 
to note that electrons in free diffusion are probably unaffected by turbu- 
lence. 

As in the case of laminar flow, it is possible to define a characteristic 
length, L s , and time, t, , for the sweeping due to the mean flow. Using 
(12), it follows that 

Lr 1 = V v /2D Te , tr 1 - V v 2 /2D Te , (20) 

where D T is the turbulent diffusion coefficient for the electrons. It can 
be shown that 

D Te = D iT \ 1 + TJTi], (21) 

where D iT is the turbulent diffusion coefficient for the ions. D T . will 
simply be written as D T for the remainder of this paper. 
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The characteristic time for electron loss by diffusion to the walls will 
also be modified by the turbulent flow. It follows that 



to. 



= D T /A, 



(22) 



III. COMPARISON WITH EXPERIMENT 

While the theoretical geometry is not duplicated by any existing ex- 
periment, some verification of the theory can be obtained from the re- 
sults of Buchsbaum and Cottingham 4 in hydrogen. Their experimental 
set-up has been sketched in Fig. 2. The results for breakdown in hydro- 
gen (H 2 ) has been plotted as a function of gas velocity in Fig. 3. The peak 
power required to produce breakdown varied from about 1.7 to 3.2 kw. 
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Fig. 2 — Geometry of the experiment of Buchsbaum and Cottingham (Ref. 



£2 

O 



o EXPERIMENTAL (REF. 4) 
THEORY 




P= 4.5 TORR 
L= 22 CM 
8= 1.1 MM 

D= 1.2CM 
S BAND 



1.0 2.0 4.0 6.0 10 20 

MEAN FLOW VELOCITY IN CM/SEC 



40 

xio- 



Fig. 3 — Comparison of the theoretical and experimental results for peak 
power as a function of flow velocity in hydrogen (experimental results of Ref. 4). 
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These results were used to calculate the effective electric field intensity 
required to produce breakdown as a function of gas velocity. In this 
calculation, the effective frequency for electron-molecular elastic colli- 
sions was taken as 2 

v c = 4.8 X HTV, sec -1 . (23) 

An approximate criteria for breakdown can be obtained from the theory 
as given by (11). (See Appendix B.) It can be shown that 

*<» = [£+£]<- (25) 

where j\ in H 2 has been measured and compared to dc data. 2 For this 
particular experiment, 4 the geometry was cylindrical with 



>-6 _ J i w i«-3 



and 



or 



t oa = 1 X 10" 6 sec, t oft = 1 X 10 sec, 



A- 



A = 0.245 cm. 

It should be noted that the length of the discharge in the direction of 
flow has been taken as 3.8 cm even though grids were not used (at y = 
and y = 3.8 cm) in the experiment. 

3.1 Laminar Flow 

When the flow of the gas is laminar, the diffusion coefficient for the 
electrons should be independent of the velocity of the flow. The magni- 
tude of the effective diffusion coefficient, D, can be calculated from the 
experimental data by noting that, for small V, 

D/A 2 » V 2 /2D for V = 500 cm sec -1 . (25) 

(This is the smallest velocity for which data is available in this experi- 
ment. See Fig. 3.) This, in turn, implies that the characteristic time for 
electron removal by diffusion to the walls is much less than the charac- 
teristic time associated with the sweeping effect at this velocity of flow. 
For the laminar case, t, is given by (17) or (18) provided the diffusion 
is ambipolar or nearly so. 

If the electron density, during the pulse off time, is reduced to a point 
where the electrons are in free diffusion, then D — » D_ , where 
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-D-P ^ 1 X 10 6 cm 2 sec -1 torr. 

As the electron density increases, D decreases. The resulting effective 
diffusion coefficient, D, , has been calculated by Allis and Rose. 7 Using 
their results and the geometry of this experiment, it follows that 

D, S 0.2ZL for n = 10 6 elec/cm 3 , 

and 

D, ^ 0.04 D_ for n = 10 8 elec/cm 3 . 

Returning to the question posed by (44) of Appendix B, we note that 
the ratio 
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n small, 


A 2 
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for 


n = 10 elec/cm , 
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for 


n — 10 elec/cm . 



Consequently, diffusion can only be neglected during the pulse on time 
when the electron density is larger than 10 6 elec/cm 3 . This would appear 
to be the case during the l-/xsec, 3-gc pulse used in this experiment. 

Since the incident microwave pulse was relatively flat, (see Fig. 2 of 
Ref. 4), the peak power can be used to calculate the effective electric 
field strength and, in turn, the ionization frequency. For a flow velocity 
of 500 cm sec -1 , the peak power required to produce breakdown at the 
end of one microsecond was 1.7 X 10 watts. This power level corresponds 
to an effective field strength of about 200 volts per cm or an ionization 
frequency of 1.5 X 10 radians per sec. Using (24), it follows that 

D/A 2 ay X 10 4 sec -1 . (26) 

This corresponds to an effective diffusion coefficient for the electrons of 
920 cm 2 sec" 1 or 

D,p^4X 10 3 cm 2 sec -1 - torr. (27) 

Using this value for D, it follows that (25) is satisfied at 500 cm sec -1 , 
since 

7 2 /2D a* 120 sec -1 . 

The magnitude of the effective diffusion coefficient for the electrons 
can be compared to the ambipolar diffusion coefficient for hydrogen as 
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measured by Persson and Brown 8 and to that calculated from the mobil- 
ity measurements of Rose. 9 The results of Persson and Brown, derived 
from afterglow data yield an ambipolar diffusion coefficient of 

D a p = 700 ± 50 cm 2 sec -1 - torr. (28) 

This is a factor of six less than the effective diffusion coefficient as given 
by (27 ) . The difference could be due to the fact that the diffusion in the 
afterglow was not ambipolar for all t, or to the presence of higher order 
modes, or to an elevated temperature for the electrons. Since the collision 
frequency was on the order of 10 9 radians per sec, it would appear that 
the plasma was rapidly reduced to thermal equilibrium during the pulse 
off time of one millisecond. It would appear reasonable to conclude 
that either or both of the first two possibilities are likely. 

The mobility data of Rose 9 yields a value for D a p which is less then 
that measured by Persson and Brown in the afterglow. 

At higher flow velocities, the effect of the flow can no longer be neg- 
lected. Equation (24) can be used together with (23, 26, 27) and Fig. 3 
of Ref . 2 to calculate the peak power required to produce breakdown in 
one microsecond. The results have been plotted in Fig. 3 for flow veloci- 
ties of less than 8000 cm/sec, together with the experimental results of 
Buchsbaum and Cottingham. 

While the agreement between theory and experiment is excellent, it 
should be noted that (24) is based on the condition 

D' V 2 
Vi ^A? + 2D" pulseon - 

From the results, it is evident that the condition is reasonable provided 
the electron density is such that the diffusion of electrons is never free, 
or provided any free diffusion mode is restricted to a negligible part of 
the pulse on time. 

IV. CONCLUSIONS AND RECOMMENDATIONS 

The conservation equation for the electron density in a laminar or 
turbulent flow (2) has been solved for the parallel plane geometry of 
Fig. 1 (11). Extension to other geometries is straightforward. 

For cases where 

v - vt - p a » D/A 2 , (29) 

it was shown that pulsed breakdown in a laminar or turbulent gas may 
be controlled by the rather simple criterion, 
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IT 2 D) 
J-XU- \2D + X>j to!! ' (30) 

where V is the mean velocity of the flow. 

The power required to produce breakdown in a turbulent gas may, 
or may not, be a function of the mean velocity of the flow. A velocity 
independent result is obtained when 

Dr/A 2 » V 2 /2D T (31) 

provided D T ~ VI, , where l e is the mixing length for the plasma, and 
provided l e is inversely proportional to the Reynold's number of the flow. 
In order to test the present theory further, future experiments should 
be designed such that 

Vi » D/A for all t. 

It would be of interest to measure the electron density in the discharge 
at a number of positions in the direction of flow. This result could be used 
to check the predicted theoretical variation as well as insure the presence 
or absence of turbulence. 
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APPENDIX A 

Equation (5) can be solved by the method of separation of variables. 
Let 

n(y,z,t) = Y(y)-Z(z)-T(t). (32) 

It follows that 

IdT ? IdY , v n /ld 2 F, ld 2 Z\ ,__. 



1 dT 1, . , 1 (d 2 Y V v dY\ . ld 2 Z /nA . 

DTM = D (pi - V ° ) +Y\W-D^yj+ZdS- (U) 

Taking 

1 d'Z 2 

ZM = ~ a 
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and using the appropriate boundary condition yields 



. m-K 



where 



Taking 



yields 



Z(z) = A m sm'-^z (35) 
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r \ 
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Y\dy> D dyf ~ P ' 
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with the solution Y(y) = e +uy , where 



w 2 _ V* u + fa m Q 



Hence, 



wm- 



«-^±^/C„M -■«.* 



If the solution is to satisfy the boundary conditions at y = and ?/ = I, 
it follows that Y(y) must be periodic in?/; hence, 4/3 2 > (7„/Z)) 2 . Con- 
sequently, 



U = 2D ±J 



with the solution 

Y(y) =B n e rv » ,2D) »sm 7 ^y (36) 

where 

Finally, the differential equation in i, reduces to 



ldT , v 



»{(?)'*(¥)'-©•}■ 
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with the solution 

T(t) -«p {(*-*) -d\(^j) r ■(■;■ ) I- (,7; 



*w+<m- <»" 



It follows that the solution of (5) satisfying the boundary conditions 
given in (4) is 



n(x,y,z,l) = €"' 2- 2-. A. mn sin — y sin — z 

m n I U 

(38) 

It is apparent at the offset that the usual method for evaluating the 
constant A mn cannot be used for the general case since 

sin -y- v/ «m . 2/ e " w d?/ (39) 

does not vanish for m 7^ n. In view of this, it is apparent that the usual 
modes do not exist. 

Fortunately, a separation into modes is still possible, provided the 
electron density, at t = 0, is given by 

f(x,y,z)t lv >> liD)v . (40) 

This would appear to be reasonable for the steady state, pulse on pulse 
off, conditions of the usual pulsed experiment. 

Let us assume that the density at t = 0, (end of the off cycle) is 
given by 

n(x,y,z,0) = m X e (v »'*»* (41) 

For this case, 

4 1 

A mn = n -. — [1 — cos nir][l — cos mir] (42) 

7T 2 mn 

from which it is evident that A „„. = for m or n even. It is of interest 
that the l 8t "mode" for this case is 

/ .\ 16 . 7T . 7T (v u l2D)y 

n{x,y,z,t) = n — sin -. y sin - z t " 

t- I a , , 

(43) 

In practice, the ratio V U /2D at any velocity will depend on the particular 
gas used. 
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For example, in the experiment of Ref. 1, 

(V/2D)y < 11 for ^ y g 3.5 cm 

in hydrogen, for V < 7 X 10 3 cm/sec, the laminar region. In the turbu- 
lent region, the ratio decreases since D increases. 

APPENDIX B 

The direct application of the theory to any particular experiment may 
be complicated by the variation of the diffusion coefficient with electron 
density and, in turn, time. The interested reader is referred to a recent 
paper by Buchsbaum and Cottingham 2 for a discussion of this problem. 
Following their example, the effects of electron attachment in H 2 will be 
neglected. 

Let the flow conditions be such that the initial density, n , is sufficient 
to insure ambipolar diffusion or near ambipolar diffusion of the electrons 
for all or nearly all time. This assumption will be justified for the par- 
ticular experiment. 

During the pulse on time the appropriate equation for the electron 
density is given by (43) with D = D' (hot electrons) and v a = 0. Let 
D' and V be such that 

„. » tL + * for all ^ t ^ t on . (44) 

A JiL) 

This implies that 

vi >5> tD~ X and t,'" 1 

for all flow velocities during the pulse on time. It follows that (43) for 
the first "mode" reduces to 

n(x,y,z,t) ^ n - — sin f 7/ sin - z « " e . (,45J 



7T 



I U " d 



This result can be applied directly to the experiment of Buchsbaum and 
Cottingham. 4 

In this experiment, electron density was monitored indirectly by ob- 
serving the intensity of the power reflected from the discharge region. 
Since the electron density and, in turn, the reflection coefficient was a 
function of position within the discharge it is evident that the experi- 
mental results can be related to (45) by choosing (x,y,z) such that the 
reflection coefficient is a maximum. Hence, 

n(t) ^ n * exp ( Vi t), for ^ / ^ t on . (46) 
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The appropriate equation for the pulse off time may be derived from 
(43 ) in a similar manner. It can be shown that 

n(t) = n mnx * exp - jp + ^| t for ^ t ;g t o!f . (47) 

Equations (46) and (47) can be solved to yield the approximate 
criteria for breakdown. 

"'™ = ! !+2dK" (48) 



LIST OF SYMBOLS 

d = distance between the parallel planes in the z direction. 
I = distance between the grids in the direction of flow. 
l e = effective mixing length for the electrons in a turbulent flow. 
li = mean free path for the ions. 
n = electron density in a laminar flow or the ensemble average of the 

electron density (at a given time) in a turbulent flow. 
t D = characteristic time for electron removal by diffusion to the walls. 
t a = characteristic time for electron removal by the flow of the gas. 
v = average velocity of the electrons or ions resulting from ambipolar 
diffusion and the mean flow. 
(vi) = mean thermal velocity of the ions. 
D a = coefficient of ambipolar diffusion for the electrons in a laminar 

flow. 
Di = diffusion coefficient for the ions in a laminar flow. 
Di T = diffusion coefficient for the ions in a turbulent flow. 
D T = diffusion coefficient for the electrons in a turbulent flow. 
L.s = characteristic "sweeping" length for electron removal by the 

flow of the gas. 
R e = Reynold's number of the flow. 
T t = temperature of the electrons. 
Ti = temperature of the ions. 

V = velocity of the gas in a laminar flow or the ensemble average of 
the gas velocity (at a given time) in a turbulent flow. Experi- 
mental conditions usually permit a time average in this case. 
a = recombination coefficient for the electrons. 
A„„, = effective diffusion distance for the ??m mode. 
v a = attachment frequency for the electrons to H + and H2 + . 
v c = elastic collison frequency for the electrons. 
v Ci = elastic collision frequency for the ions. 
Vi = ionization frequency for the electrons. 
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